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Abstract 

In this paper, the concepts and the direct theorems of stability in 
the sense of Liapunov, within the framework of Birkhoffian dynamical 
systems on manifolds, are considered. The Liapunov-type functions are 
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1 Introduction 

During the last years, a far reaching generahzation of the Hamiltonian frame- 
work has been developed in a series of papers. This generalization, which is 
based on the geometric notion of generalized Dirac structure (see Courant [S] 
and Dorfman [6J), gives rise to implicit Hamiltonian systems (see, for example, 
the papers by Maschke and van der Schaft [12], [14]). Applications to non- 
holonomic systems and electrical circuits (see Bloch and Crouch [2], Maschke 
and van der Schaft [H]) illustrate this theory. Recently, the notion of implicit 
Lagrangian system has been developed by Yoshimura and Marsden |16j . Non- 
holonomic mechanical systems and degenerate Lagrangian systems such as LC 
circuits can be systematically formulated in the implicit Lagrangian context in 
which Dirac structures are also used. 

An alternative approach to the study of dynamical systems is the Birkhoffian 
formalism. This is a global formalism of implicit systems of second order ordi- 
nary differential equations on a manifold. It applies to a wide class of systems, 
among them, nonholonomic systems, degenerate systems as well as dissipative 
systems. Kobayashi and Oliva developed in ^ the framework of Birkhoffian 
dynamical systems on manifolds, following Birkhoff 's ideas presented locally in 
[Tj. The space of configurations is a smooth m-dimensional differentiable con- 
nected manifold and the covariant character of the Birkhoff generalized forces is 
obtained by defining the notion of elementary work, called Birkhoffian, a special 
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PfafRan form defined on the 2-jets manifold. The dynamical system associated 
to this Pfaffian form is a subset of the 2-jets manifold which defines an implicit 
second order ordinary differential system. The notion of Birkhoffian allows the 
introduction of the intrinsic concepts of reciprocity, regularity, affine structure 
in the accelerations, conservativeness [5], dissipativeness [3- 

The Birkhoffian formalism in the context of electrical circuits was discussed 
by lonescu and Scheurle [8] for the case of LC circuits, and lonescu [7j for 
the case of RLC circuits. An LC/RLC circuit, with no assumptions placed on 
its topology, will be described by a family of Birkhoffian systems, parameter- 
ized by a finite number of real constants which correspond to initial values of 
certain state variables of the circuit. It is shown that the Birkhoffian system 
associated to an LC circuit is conservative. Under certain assumptions on the 
voltage-current characteristic for resistors, it is shown that a Birkhoffian system 
associated to an RLC circuit is dissipative. For LC/RLC networks which con- 
tain a number of loops formed only from capacitors, the Birkhoffian associated 
is never regular. A procedure to reduce the original configuration space to a 
lower dimensional one, thereby regularizing the Birkhoffian, is presented as well. 

For RLC electrical networks, Brayton and Moser [3] proved under a special 
hypothesis, that there exists a mixed potential function which can be used to 
put the system of differential equations describing the dynamics of such a net- 
work, into a special form (see §4 in [3]). The hypothesis they made is that the 
currents through the inductors and the voltages across the capacitors determine 
all currents and voltages in the circuit via Kirchhoff 's law. The mixed potential 
function is constructed explicitly only for the networks whose graph possesses a 
tree containing all the capacitor branches and none of the inductive branches, 
that is, the network does not contain any loops of capacitors or cutsets of in- 
ductors, each resistor tree branch corresponds to a current-controlled resistor, 
each resistor co-tree branch corresponds to a voltage-controlled resistor (see §13 
in [3]). Making different assumptions on the type of admissible nonlinearities in 
the circuit, this mixed potential function is used in [3] to construct Liapunov- 
type functions to prove stability. 

Smale also develops the differential equations for nonlinear RLC electrical 
circuits and illustrates these equations through a series of examples. He builds 
on the work of Brayton and Moser [3] but he is able to treat more general equa- 
tions. A large part of the paper illustrates these equations by means of examples 
and discusses stability properties of the examples. 

The paper at hand is organized as follows: at the beginning of Section 2 
we present the basics of Birkhoffian systems, from the viewpoint of differential 
geometry using the formalism of jets. Then, we introduce the concepts and 
the direct theorems of stability in the sense of Liapunov, within the framework 
of Birkhoffian systems. In Section 3 we consider, in turn, the linear and the 
nonlinear LC networks, as well as the linear and the nonlinear RLC networks. 
For each of them we construct Liapunov-typc functions to prove stability or 
asymptotic stability under certain conditions. Finally, we discuss an example 
in Section 4. 
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2 Liapunov's direct method for BirkhofRan sys- 
tems 



In order to present the ideas in a coordinate free fashion, we consider the for- 
mahsm of 2-jets. Let M be a m-dimensional differentiable connected manifold. 
We consider the tangent bundles (TM, ttm, M) and {TTM, tttm, TM). 
The 2-jet bundle (J'^iM), ttj, TM) is defined by 

J2(M) := {z e TTM / T^m{z) - ^tm{z)} (2.1) 

where {Tttm)v ■ TyTM — > Ttt^j-^^-jM is the tangent map and 

TTj := 7rTM|j2(M) = TTTMlj-^iM) (2.2) 

A local system of coordinates (q) — {q-')j=i,....m on M induces natural local co- 
ordinates on J'^{M), denoted by {q,q,q) = {q^ , q^ , q^)j=i m (see for example 

A Birkhoffian corresponding to the configuration manifold M is a smooth 
1-form uj on ,P{M) such that, for any x g M, we have 

i> = (2.3) 

where Lx ■ (3~^{x) — > J'^{M) is the embedding of the submanifold (3~^{x) into 
J^(M), (3 = TTM ° TTj. From this definition it follows that, in the natural local 
coordinate system {q, q, q) of J^{M), a Birkhoffian oj is given by 

m 

u; = ^Q,(g, 9, g)dg^ (2.4) 
i=i 

with certain functions Qj : J'^{M) R. The pair (M, uj) is said to be a 
Birkhoffian system (see [5]). 

The differential system associated to a Birkhoffian lu is the set (maybe empty) 
D{uj), given by 

D{uj) -.^ {z £ J^{M) \lu{z) ^0} (2.5) 

The manifold M is the space of configurations of D{u)), and D(uj) is said to have 
m 'degrees of freedom'. The Qi are the 'generalized external forces' associated 
to the local coordinate system. In the natural local coordinate system, D{lli) is 
characterized by the following implicit system of second order ODE's 

Oj (g, g) = for aU j = (2.6) 

The Birkhoffian formalism is a global formalism for the dynamics of implicit 
systems of second order differential equations on a manifold. 

A cross section X of the affine bundle (J^(M), ttj, TM), that is, a smooth 
function X : TM — > J'^{M) such that TrjoX^id, can be identified with a special 
vector field on TM, namely, the second order vector field Y on TM, that is, a 
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smooth function Y : TM -> TTM such that tttm o Y=id and Tttm o Y =id. 

Using the canonical embedding i : J^{M) — > TTM, we write Y — i o X . 

In natural local coordinates, a second order vector field can be represented as 



d 



d 



(2.7) 



A Birkhoffian vector field associated to a Birkhoffian lu of M (see 9 ) is a smooth 
second order vector field on TM, Y ^ioX, with X : TM J^M), such that 
ImX C D{ijj), that is, 

X*uj = (2.8) 

In the natural local coordinate system, a Birkhoffian vector field is given by the 
expression (|2.7p . such that £lj{q, q, q{q, q)) — 0. 
A Birkhoffian lu is regular if and only if 



det 



50, 



dqi 



{q, q, q) 



^0 



(2.9) 



- J,?— l,...,m 

for all (g, q, q), and for each (q, g), there exists {q,q,q) e J'^{M) such that 
9, 9) = 0, i = 1, ...,TO. 
If u; is a regular Birkhoffian corresponding to the configuration manifold M, 
then, the principle of determinism is satisfied, that is, there exists an unique 
Birkhoffian vector field Y ^ i o X associated to lu such that ImX — D{lu) (see 
i). 

A Birkhoffian uj of M is called conservative (see [9]) if and only if there exists 
a smooth function E^, : TM R such that 



{X*lu)Y = dE^{Y) 



(2.10) 



for all second order vector fields Y ^ i o X , which is equivalent, in the natural 
local coordinate system, to the identity 



dE^ 



dqi 



(2.11) 



If UJ is conservative and F is a Birkhoffian vector field, then (j2.10|) becomes 

dE^{Y)^Q (2.12) 

This means that E^ is constant along the trajectories of Y . 

A Birkhoffian lu of the configuration space M is called dissipative (see [7]) if 
and only if there exists a smooth function Eq^ : TM R such that 



{X*Lu)Y ^dEojY)+D{Y) 



(2.13) 



for all second order vector fields Y = ioX on TM, D being a dissipative 1-form 
on TM, that is, D ^ X^jli Dj{q,q)dq^ and 



(2.14) 
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Equation (j2.13p is equivalent, in a local coordinate system, to the identity 



7n m 

In view of p.l4p . we obtain from (|2.13p . 

{X*Lu)Y >dEo^iY) (2.16) 

for all second order vector fields Y ^ i o X . That is equivalent, in local coordi- 
nates, to the dissipation inequality 

m m 
^0, (<Z q, iiW > 

If w is a dissipative Birkhoffian and Y is the Birkhoffian vector field, then (|2.16p 
becomes 

dEQ^(Y)<0 (2.18) 

This means that Eq^ is nonincreasing along the trajectories of Y. 

Let us introduce now the concepts of stability for a Birkhoffian system. 
The equilibrium points of the system, that is, the points in which the system 
can remain permanently at rest, are to be found as the solutions of the system 

0,((?,O,O) = 0, j = l,...,m (2.19) 

Let us denote an equilibrium point by (ge, 0) G C TM, and an initial state of 
the system by (go, 90), with q{0) go, 4(0) = qo- 

For regular Birkhoffians, we can define the equilibrium points using the notion 
of Birkhoffian vector field, that is, a point {qe, 0) is an equilibrium point of the 
Birkhoffian vector field Y if and only if 

Y{q,,0) = (2.20) 

An equilibrium point (qg , 0) is said to be stable (or Liapunov stable) if for 
every open neighborhood of (qe, 0), there is a neighborhood fii C f2 such that 
a motion {q{t),q{t)) starting at (go,9o) G ^^i, remains in fl. If in addition, fli 
can be chosen such that, for any (go, 9o) G ^^i, {lit), q{t)) converges to (qe, 0) as 
t — > 00, then {qe, 0) is said to be asymptotically stable. 

In the memoir '11', Liapunov presents geometric theorems, generally referred 
to as the direct method of Liapunov (see, for example, [TO]), for deciding the 
stability or instability of an equilibrium point of a differential equation. 
In what follows we consider Liapunov 's direct method for Birkhoffian 
systems. This is based on finding a function F £ {TAI, R) such that 

(j) F(ge,0) = 

(il) Viq, q)>0 for {q, q) ^ {q^, 0) in O (2.21) 
(Hi) dV{Y) < for all second order vector fields Y defined on 



dEo^ 



dqi 



-q' +D,{q,q)q^ 



(2.15) 



dEp^ 
dqi 



dqi 



(2.17) 
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with ri an open neighborhood of ((7e,0). The function V is cahed Liapunov 
function. 

One can prove the following theorems (completely analogous to the theorems 
proved in lOj for a Liapunov function defined on C/ C M): 

Stability theorem. If there exists in a neighborhood fl of {qe,0) a Lia- 
punov function V{q,q), then (ge,0) is stable. 

Asymptotic stability theorem. If there exists in a neighborhood of 
(ge,0) a Liapunov function V{q,q), such that dV{Y) < for all second order 
vector fields Y defined on il, then {qe,0) is asymptotically stable. 

From the condition (ii) in (|2.2ip we get that there exists Cq > such that 
the level curve 

{iq,q)en,V{q,q)=c} (2.22) 

is a closed curve for every constant < c < cq. Sketching in the m-plane {q,q) 
these level curves of the function V, we obtain surfaces like "ellipsoids" centered 
at the equilibrium point. 

If dV = 0, then the equilibrium point (<Ze,0) is a center and the motion of the 
system is periodic. 

If dV < 0, then each trajectory keeps moving to lower c and hence penetrates 
smaller and smaller "ellipsoids" as < — > cw. Thus, the equilibrium point is asymp- 
totically stable. This exclude the existence of periodic motions of the system. 

3 Stability of the equilibrium points of LC and 
RLC networks 

A simple electrical circuit provides us with an oriented connected graph. The 
graph will be assumed to be planar. Let b be the total number of branches in 
the graph, n be one less than the number of nodes and m be the cardinality of 
a selection of loops that cover the whole graph. By Euler's polyhedron formula, 
b = m + n. We choose a reference node and a current direction in each Z-branch 
of the graph, I — We also consider a covering of the graph with m 

loops, and a current direction in each j-loop, j — \,...,m. We assume that 
the associated graph has at least one loop, meaning that m > 0. An oriented 
connected graph can be described by matrices which contain only 0, ±1, these 
are: the incidence matrix B £ 9Jlb„(R), rank(i3) = n, and the loop matrix 
A G 9Jlbm(R-), rank(A) — m. For the fundamentals of electrical circuit theory, 
see, for example, [4]. 

Let us now consider an RLC electrical circuit consisting of r resistors, k 
inductors and p capacitors, such that to each branch of the associated graph 
there corresponds just one electrical device, that is, b — r + k + p. For LC 
electrical circuits r — 0. Using the matrices A and B, Kirchhoff's current law 
and Kirchhoff's voltage law can be expressed by the equations 

B^i^O (KCL), A^v^O (KVL) (3.1) 
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where i — (i[r]7i(a))ia) G ^ x x R-^ — R-'' is the current vector and v = 
{^\r\-,V{a)i'Va) e R*" X R*^ X RP ~ R*" is the voltage drop vector. Tellegen's 
theorem estabUshes a relation between the matrices and : the kernel of 
the matrix is orthogonal to the kernel of the matrix A"^(see, for example, 
page 5 of [5]). 

We consider the voltage-current laws for nonlinear devices given by 

■yp] = i?r(l[r]), V(^a) = La{l(a))—^, Va=Ca{Qa), (3.2) 

Rr,La,Ca ■ R — '' R\{0} being smooth functions, Q^, denote the charges of 
the capacitors, with Iq = If the capacitors and the inductors are linear 

then the relations above become, respectively, 

v[r] = Rri[r], f(a) = ^^—j— Va = —, (3.3) 

where Rp 7^ 0, 7^ and 7^ are distinct constants. 
Summing up, the equations governing the network are 

/ ^r(i[r]) \ 
= 0, ia(i(.))%i =0. (3.4) 

Using the first set of equations (13. 4|) . one defines (see [7j, [8]) a family of m- 
diniensional affine- linear configuration spaces C R*", parameterized by a 
constant vector c in R" which corresponds to initial values of certain state 
variables of the circuit. Since the matrix B is constant, integrating the first set 
of equations (j3.4p . one gets B'^x — c, with i = i, c a constant vector in R". 
Thus, one defines 

M, := {a; £ Il''\B^x = c} (3.5) 

Its dimension is m = b — n, because rank(i3) — n. Local coordinates on are 
denoted by q = {q^,..,q^"). Solving the system in (|3.5p . one expresses any of 
the a;- variables in terms of q-s, namely, 

x=Afq + IC (3.6) 
where J\f = I A/"" I is a matrix of constants and K. = 





/C" I a constant vector in R*" 

By Tellegen's theorem and a fundamental theorem of linear algebra, one obtains 
that KeriA'^) = Ker{Af^) (see [7] , [5]). 

A Birkhoffian uj^ on the configuration space arises from a linear combination 
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of the second set of equations (|3.4p . by replacing the matrix with the matrix 
of constants Af'^ . 

I) For a hnear LC network (r = 0) we have the foUowing expression of the 
Birkhoffian (see [8j) 

km b m J\faj\fa 

Q,(<z,9,g)-EEL--^°-^°«^+ E E ^^^'^^ + (3.7) 

1-1 7 11-1 ^a—k 

with const £ R™ a constant vector. 

A hnear LC network is conservative (see [5]). The function : TMc — > R 
satisfying (|2.1ip has the following expression 



a—1 1 a— K+l 1 

m 

-E(const)j(7^ (3.8) 



In what follows we assume that 

k 



det 



ELaA/?A/^" 



a=l 



^ 0, det 

l,...,m 



.a— fc+l 



^ (3.9) 

j,i—l....,m 



that is, the network does not contain loops formed only by capacitors and re- 
spectively, loops formed only by inductors (see [8]). If the network contains 
capacitor loops and inductor loops, we will reduce first the configuration space 
to a lower dimensional configuration space. On the reduced configuration space 
the corresponding Birkhoffian is still conservative (see [5]) and the correspond- 
ing determinants (|3.9p will be different from zero. The inductor loops can be 
considered as some conserved quantities of the network. 

Theorem 1. Let (qe, 0) be an equilibrium point of a linear LC network with 
the Birkhoffian components given by |5'. 7] j. Then qe satisfies the system 

E E^^ + (const)j = 0, j = l,...,m (3.10) 

J , 1 ■ 1 ^a~k 

To each const which is related to the initial data for the considered network, we 
get a unique equilibrium point. If 

La > 0, Va = 1, ...,A:, > 0, Va = 1, (3.11) 

the equilibrium point is a stable center, and the motion of the system is peri- 
odic. 

Indeed, the equilibrium points of a linear LC network are obtained as solu- 
tions of the system O.j{q,0,0) = 0, j — 1,...,to, where 0.j{q,q,q) is given by 
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p.7p . Thus, we get that Qe has to fulfill the system (|3.10p . Under the second 
condition in (|3.9[) . this system has for each const S R™ a unique solution. 
The stability of this equilibrium point is obtained by the Stability Theorem 
presented in section 2. We define a Liapunov function V G C^{TMc, R) by 

= 2 E E ^aM^'KQ'<i' + 2 E E - - 

a— 1 1 a—k-\-l 

(3.12) 

where ge satisfies the system p.lOp . Indeed, this function satisfies the condi- 
tions (|2?2T|) . Taking into account (l3TT|) . the matrices (I]q=i LaA^fTV^") and 

(^''g.^f^^i positive definite. Thus, the condition (m) in (|2.2ip is 

fulfilled. The first determinant in (|3.9p being different from zero implies that 
the corresponding Birkhoffian is regular. Therefore, along the trajectories of 
the unique (principle of determinism) Birkhoffian vector field, the function 
defined in (|3.8p . satisfies ()2.12|) . Thus, the function (I3.12p satisfies the condition 
(iii) in ()2.2ip . In this case, sketching in the m-plane {q,q) the level curves of 
the function (|3.12p . we obtain ellipsoids centered at the equilibrium point. The 
equilibrium point is a center and the motion of the system is periodic. □ 

II) For a nonlinear LC network we have the following expression of the 
Birkhoffian (see [8]) 

k / m \ m b / m \ 

0,(g, q, q) = Y.KLa E-^??' E-^"9' + E -^/c-.., Y.^ri' + 

a^l ) z^l a^/c+l / 

m / k \ h 

= E E-^/^"^- ^i) + E ^tCo.-k {q) (3.13) 

i=l \a=l ) Q=fc+1 

A nonlinear LC network is conservative (see [8]). In this case, the function 
: TM, ^ R is given by 



E^iq,q)=£{q) + iBiq) (3.14) 



with 



k m m r r 

^i<i) = EE E (-1^7 [Lii-'\<i)m 



a— 1 / — I Zi<...<Z;— 1 



+ {l-\)Lt''\q^Ml...Mld<tK..dcl^ 

<J lib lit n n 

E E E (-1)'^7 \&t'^mt...M-dq^K..dq\i.\^) 



a=fc+l i=l ii<...<i, = l 
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where we denoted the derivatives C^'lj. := 



In what follows we assume that 
fc 



det 



.0=1 



^0 (3.16) 



that is, the network does not contain capacitor loops. In the case the network 
contains capacitor loops, we first reduce the configuration space to a lower di- 
mensional one. On the reduced configuration space the corresponding Birkhof- 
fian is still conservative (see |S ) and the corresponding determinant above will 
be different from zero. 

Theorem 2. Let {qe,0) be an equilibrium point of a nonlinear LC network 
with the Birkhoffian components given by 113.13]} . Then satisfies the system 

b / m \ 

5] A/-/C„_J^AArq'+/C" =0, j-l,...,m (3.17) 

a=k+l \l=l / 

A nonlinear LC networks can have several equilibrium points. If 

La(0)>0, Va = l,...,fc, C^(ge) >0, Va = l,...,p (3.18) 

then the equilibrium points are locally stable centers. 

Indeed, the equilibrium points of a nonlinear LC network are obtained as 
solutions of the system 0j(g,O, 0) — 0, j ^ 1, ...,m, where £}j{q,q,q) is given 
by (|3.13[) . Thus, we see that q^ has to fulfill the system (|3.17[) . 
The local stability of the equilibrium points is obtained by the Stability Theorem 
presented in section 2. We define a Liapunov function V £ C^{TMc,IV) by 

V{q,q)^Em<i)-EUq,,0) (3.19) 

with given by (|3.14|) and qe satisfying the system (|3.17|) . 
Let us now evaluate the Hessian matrix of the function V in (|3.19|) at an equi- 
librium point {qe,0). We get 







nv{q.,0)-\ (3.20) 

For the Birkhoffian ([XT^ . the function in (|XTI)) satisfies the identity (piTj) 
(see [H]), that is, 



^ - EE^^"(^)W^' (3-21) 

a=l 1 = 1 
b 

J2 C^-k{q)K'' (3.22) 

a=k+l 



dq 

d^<Biq) 
dq' 
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Therefore, we get 

d'^eiq) 



J2 E Kim.^KK'i' + Laim^M^ 

a=l L/=l 
b 

a=k+l 



From dnSSl), (ESI), the matrix ([X^ writes as 



Hy(ge,0) = 



E'=iia(0)A/;°-AA' 



(3.23) 
(3.24) 

(3.25) 



and 



ELfe+iC^4-.(9e)A/;"A/;" 

In view of the conditions (|3.18p . the matrices |^X]a=i (0)-A/^''A/'° 

(X^Q^fc+i ^Q-fe('?e)-A/^"-A/'"^ are positive definite. Therefore, the Hessian ma- 
trix (|3.25p is positive definite. The centers of the level curves of the function 
p.l9p have the coordinates ((7e,0), where q^ satisfies the system p.l7p . Thus, 
in a neighborhood of an equilibrium point, the condition {ii) in (j2.2ip is fulfilled 
by the function V in (|3.19p . The determinant (|3.16p being different from zero 
implies that the corresponding Birkhoffian is regular. Therefore, along the tra- 
jectories of the unique (principle of determinism) Birkhoffian vector field, the 
function Euj satisfies p.l2p . Thus, the function (|3.19p satisfies the condition 
{Hi) in (|2.2ip . By the Stability Theorem, the equilibrium points are locally 
stable centers. □ 



III) For a linear RLC network we have the following expression of the 
Birkhoffian (see [7]) 



r+k m 



£l,iq,q,q) = E E ^a--^'^-^"?' + E E «r-<-^f 9* 



a—r-\-l i—1 
b 



r=l i=l 



E E^'^' + (const), 



(3.26) 



a=r+k+l i=l 

with const £ R™ a constant vector. 
A linear RLC network with 

Rr > 0, r = 1, .. 

is dissipative (see [7]). The function Eq^ : TAI^ 
satisfying (I2.15p . are given by 

r-\-k m 



(3.27) 

R and the dissipative 1-form 



^ I ] rL III. -| U IIL \ I 

EoMA) = ^E H^^-r^mi'i'^ o E Ett^ 



a— r+1 i,*— 1 
in 

E (const) J- 



^ fi — ?^ — k 

a.—r-\-k+l 



(3.28) 
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E E^rA/J-^^W (3.29) 

j,j=ir=i 

In what follows we assume that 



det 



r+k 



a—r-\-l 



^ 0, det 

-I l,...,m 



^0 (3.30) 

j.i— l,...,m 



that is, the network does not contain capacitor loops and inductor loops, re- 
spectively. If the network contains capacitor loops and inductor loops, we will 
reduce first the configuration space to a lower dimensional configuration space. 
On the reduced configuration space the corresponding Birkhoflian is still dissi- 
pative (see [7|) and the corresponding determinants above will be different from 
zero. 

Theorem 3. Let {qe,0) be an equilibrium point of a linear RLC network 
with the Birkhojfian given by i3.26]) . Then satisfies the system 

b ™ A/" "TV" 
E E^ + (const), = 0, i = l,...,m (3.31) 

, , , T - T Ca — r — k 
a—r-\-k+l i—i 

To each const whieh is related to the initial data for the considered network, we 
get a unique equilibrium point. If 

La > 0, Va = l,...,fc, > 0, Va = 1, (3.32) 

the equilibrium point is asymptotically stable. 

Indeed, the equilibrium points of a linear RLC network are obtained as solutions 
of the system £lj{q, 0,0) — 0, j — 1, ...,to, where £lj{q,q,q) is given by ()3.26|) . 
Thus, we see that q^ has to fulfill the system (I3.31|) . Under the second condition 
in p.30p . this system has for each const G R™ a unique solution. 
The asymptotic stability of this equilibrium point is obtained by the Asymptotic 
Stability Theorem presented in section 2. We define a Liapunov function V € 
C\TM,,R) by 

^ r-\-k m 

V{q,q) - E^^iqA) - E^^{q,,0) = - E "^--rKKi' 

a— r+1 J, 1—1 

1 ™ TV" TV" 

+ 9 E E 77^^(9' (3-33) 

a—r+k+1 J, 2—1 

where qe satisfies the system p.3ip . Indeed, this function satisfies the conditions 
(H^. Taking into account ([O^ . the matrices (Eltr+i La-,TV°TV^") and 

fSa=r+fe+i ) positive definite. Thus, the condition {ii) in (|2.2ip 

is fulfilled. The first determinant in (|3.30p being different from zero implies 
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that the corresponding Birkhofhan is regular. Therefore, along the trajectories 
of the unique (principle of determinism) Birkhofhan vector field, the function 
Eq^ satisfies (j2.18p . Thus, the function (|3.33p also satisfies (|2.18p . In this case, 
sketching the level curves of the function (|3.33p in the m-plane (g, g), we ob- 
tain ellipsoids centered at the equilibrium point. From the Asymptotic Stability 
Theorem we conclude that the equilibrium point is asymptotically stable. This 
excludes the existence of periodic motions of the system. □ 

IV) For a nonlinear RLC network we have the following expression of the 
Birkhoffian (see [7|) 

Q, (g, g, q) = ^ ^J^-- E"^"?' E-^"^' + 

r / in \ h / m \ 

r^l / a^r+A;+l / 

m r-i-k r b 

= E E A/7-^"ia-.(<z)9* + E-<^r(g)+ E ^r^a-r-k{q) 
i=l a=r+l r=l a=r+k+l 

(3.34) 

In order to obtain a dissipative Birkhoffian (see [7]), we assume that, for all 

xRr{x)>0, Vr = l,...,r (3.35) 

that is, for each nonlinear resistor, the graph of the function Rr lies in the first 
and third quadrants. The function Eq^ : TM^ — > R and the dissipative 1-form 
satisfying (|2.15p . are given by 

Eo^{q,q)=£oiq) + €o{q) (3.36) 

with 



I T^rv III, lit I* I* 

^o(q) - E E E (-1)'^' / / [^tPiq)m 



+ il-l)L^^:^\q)\M.^^...N-t,dr-dq" 

b m 7n r r 

^oiq) = E E E (-l)'""'/ ct^^{q)K-K,dq^'-dm^) 



a=r+k+l 1=1 ii<...<ii = l 



and 



^ = EE-<^r((7)dg^' (3.38) 

3=1 r=i 
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In what follows we assume that 



det 



r-\-k 

E 

_a— r+l 



J2 La^MNfMt 



^ (3.39) 



that is, the network does not contain capacitor loops. In the case that the 
network contains capacitor loops, first, we first reduce the configuration space to 
a lower dimensional configuration space on which the corresponding Birkhoffian 
is still dissipative (see ^) and on which the corresponding determinant above 
will be different from zero. 

Theorem 4. Let {qe,0) be an equilibrium point of a nonlinear RLC network 
with the Birkhoffian given by {3.34% Then q^ satisfies the system 

r h / m \ 

^A/7i?r(0)+ A/'/Ca-r-fc E-^"9'+'^" =0, j = l,...,m (3.40) 

r=l a=r+k+l \l=l ) 

A nonlinear RLC network can have several equilibrium points. 

1) // 

i?r(0) = 0, Vr = l,...,r (3.41) 
La(0)>0, Va = l,...,fc, C;,(ge) >0, Va-l,...,p (3.42) 
the equilibrium points are locally asymptotically stable. 

2) // there exists F = 1, r such that -Rr(O) ^ 0, but, for all x ^ Q, 

x(i?r(a;) -i?r(0)) > 0, Vr = l,...,r (3.43) 

and the conditions {3.4^2^ are fulfilled, then the equilibrium points are locally 
asymptotically stable. 

Indeed, the equilibrium points of a nonlinear RLC network are obtained as 
solutions of the system 0j(g,O, 0) = 0, j = 1, ...,m, where 0.j{q,q,q) is given 
by p.34p . Thus, we see that qe has to fulfill the system (I3.40p . The local 
asymptotic stability of the equilibrium points now follows from the Asymptotic 
Stability Theorem presented in section 2. 

First we assume condition (|3.4ip to be satisfied. Then the system (|3.40p writes 
as 

b / m \ 

AA/C„_,_J^AAr(z'+/C" =0, j = l,...,m (3.44) 

a=r+k+l \l=l / 

In order to show 1), we define a Liapunov function V £ C^{TM^,'R) by 

Viq,q)^Eo^iq,q)-Eo^iq,,0) (3.45) 

with Eq^ given by (|3.36p and qe satisfying the system p. 441) . In the neigh- 
borhood of any equilibrium point, this function satisfies the conditions (|2.21|) . 
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Taking into account the conditions (|3.42p . the Hessian matrix of the function V 
in (|3.45p . at the equihbrium point (qe,0), 







(3.46) 

is positive definite. The centers of the level curves of the function (|3.45p have the 
coordinates (qe, 0), where qe satisfies the system (|3.44l) . Thus, in a neighborhood 
of an equilibrium point the condition (ii) in (j2.2ip is fulfilled by the function 
V in (|3.45p . The determinant (|3.39p being different from zero implies that the 
corresponding Birkhoffian is regular. Therefore, along the trajectories of the 
unique (principle of determinism) Birkhoffian vector field, the function Eq^ sat- 
isfies (|2.18l) . Thus, the function (|3.45p also satisfies (|2.18p . By the Asymptotic 
Stability Theorem, the equilibrium points are locally asymptotically stable. 
We assume now that there exists F = l,...,r such that i?.r(0) ^ 0. In or- 
der to show 2), we consider instead of the function Eq^ the following function 
Eo : TM, R 



Eo„(g,g) = £o{q) + ^0(9) + J2Aff Rr{0)q' 



(3.47) 



r=i 



Soiq), 20(9) being given by (|3.37p . and instead of D the following dissipative 
1-form 



j=i r=i 



Rt (g)-i?r(0) 



dq^ 



(3. 



In view of assumption (|3.43p . the vertical 1-form (|3.48p is indeed dissipative, 
that is. 



j=i r=i 



RriY.JVfq']^RriO) 



> 



(3.49) 



One can easily check that for the function Eq^ (g, q) given by (|3.47p and the dis- 
sipative 1-form p.48p . the Birkhoffian (|3.34p is dissipative, that is, the identity 



^0j(g q, iiW 



dEo^ aEo, 

—q H 

dq^ dq3 



Dj(g,g)g-' 



(3.50) 



is fulfilled. 

We define now a Liapunov function V G C^{TM^,'R) by 
V(g, q) = Eo^ (g, q) - Eq^ {qe, 0) 



(3.51) 



If w is a dissipative Birkhoffian and Y is the Birkhoffian vector field, then (|2.18p 
becomes 

dEo„(y)<0 (3.52) 
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The function V in (|3.5ip satisfies p.52p as well. 

The centers of the level curves of the function (|3.5ip have the coordinates (ge, 0), 
where qe satisfies the system p.40p . The Hessian matrix of the function V 
in (|3.5ip has at the equilibrium point the same expression p.46p . By the 
Asymptotic Stability Theorem, the equilibrium points are locally asymptoti- 
cally stable. □ 



4 Example 

We consider an electrical circuit with the associated oriented connected graph 
as in Figure 1. 




We have r — 1, k = 2, p = 3, n = A, m = 2,b = 6. Wc choose the reference 
node to be V5 and the current directions as indicated in Figure 1. We cover the 
associated graph with the loops Ii, The branches in Figure 1 are labelled 
as follows: the first branch is the resistive branch Ri , the second and the third 
branch are the inductive branches Li, L2 and the last three branches are the 
capacitor branches Ci, C2, C3. The incidence and loop matrices, B e 9}t64(R) 
and A e 9Jt62(R)j write as 



B = 





1 


-1 








\ 













1 









1 


-1 












1 





-1 






-1 













v 








1 


-1 


/ 



A = 



I 1 
1 


1 
1 

V 



\ 


1 

-1 


1 / 



(4.1) 



One has rank(i3) = 4, rank(A) = 2. Kirchhoff's current law and Kirchhoff's 
voltage law can be expressed by the equations 



B^\ = [KCL), A^v = {KVL) 



(4.2) 



where I = (l[r],l(a),ia) G R x x and v 
is the voltage drop vector. 

We define the configuration space by 



{V[T],V(a),Va) G R X R^ X R^ 



61 dT„ 



M, := {x e IC\B 



c} 



(4.3) 
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with c a constant vector in R'*. A^f^ is an affine-linear subspace in R^, its 
dimension is 2. The system in (|4.3p writes as 



-x^ + x^ 



Cl 
C4 



(4.4) 



We denote local coordinates on by q — {q^ , q"^) . If we take, for example, 

.1 .5 „2 .6 (4 5) 



x^,q^ := x^ 



we get 



q 



Cl 

Cl - 

C3 

„2 



C2 + C3 + C4 



(4.6) 



Cl + C2 



C3 



IC 



Thus, the matrix of constants Af in p.6p is exactly the matrix A and the constant 

Cl + C2 + C3 + C4 
-C3 
Cl + C2 + C3 



J 

First we consider the case that all the electrical devices in the circuit are 
linear, they are described by the relations (|3.3p . In this case, in terms of the 
g-coordinates (j4.5p , the Birkhofhan uj^ on writes as 



V 



Riq 

1 



Cl 



-q 



1 

— 1 




Cl 


C2, 




1 


vCl 


C3 



1 

— ( 

Cl 

2 



Cl + C2 + C3 



Cl 



C2 + C3 



Cl 



(4.7) 



Let us see now how the constants are related to the initial conditions that 
may be specified for the considered network. 

The differential system associated to the Birkhoffian ()4.7p is written 



Qi(<Z,g,g) =0, £l2{q,q,q)^0 



(4.8) 



For each capacitor we are able to specify the initial charge, that is, Qi(0), Q2(0), 
Q3(0), and for each inductor the initial current, that is, i(i)(0), i(2)(0). Taking 
into account (|4.5p . the relation i = i and the relations two and three in (|4.7p . 
we have the following initial conditions for the differential system (|4.8p 



q\0) = Q2(0) 
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q\0) = Q3(0) 
q\0) = l(i)(0) 

e{0) = 1(2) (0) (4.9) 

Besides, taking into account the notations (|4.5p and the last relation in (14. 7p . 
we find 

Cl + C2 + C3 = Qi(0) - Q2(0) + Q3(0) (4.10) 

Thus, the Birkhoflian (14.71) becomes 



/ • -N -1 -1/1 1\ 1 1 2 Qi(0)-Q2(0) + Q3(0) 
0i(g,(?,g)-Ligi+Rigi+ — + — ^2 , viv ; ^2V y W3V ; 



. - 1,. + (i- + 1) ,^ - + (411) 

Cl \Ci CsJ Cl 

If the constant Ri > 0, the BirkhofSan (|4.1ip is dissipative. The function 
Eq^ : TM^ R and the dissipative 1-form satisfying (|2.15p have the expressions 

Eo.iqA) - ^^i(^^)^ + ^^2(^^)^ + ^(^^-^^)^ + ^(^T + ^(.^f + 

Qi(0)-Q2(0) + Q3(0) 1 _ Qi(0)-Q2(0) + Q3(0) 2 .4 ^2) 

Cl ^ Cl 

D = Ridq^ (4.13) 

The equilibrium point of the considered linear network is the solution of the 
system 

1 n 1 1 2 . Qi(0)-Q2(0) + Q3(0) n 
_l^i ^ n ^ 1\ 2 _ Qi(0)-Q2(0) + Q3(0) ^ (4,^4) 

Cl VCl C3y Cl 

If the constants Li, L2, Ci, C2, C3 satisfy the conditions p.32p . this equilibrium 
point is asymptotically stable. We define a Liapunov function V by 

V{q,q) = Eo^iq,q)-Eojqe,0) = ^Li{q')^ + ^L2{qy 

+ ^ - q') - ill - + ^(-Z^ - lir + ^{q' - qll^-l^) 

ZCi zC2 ^C3 

where q^ satisfies the system (|4.14p . The level curves of the function (|4.14p 
represent a set of ellipsoids surrounding the equilibrium point. Because the 
Birkhoffian (|4.1ip is dissipative, it follows that dE^^ < 0, and therefore dV < 0. 

Let us consider now the case that all the devices are nonlinear, they are 
described by the relations ()3.2p . For the coordinate system on Mc given by 
(|4.5p . the Birkhoffian becomes 

0i(g,g,g) = Li(gi)gi + i?i(gi) + C,{q' - q' + IC') + C^iq') 
02(g,<7,9) = L2{q')f - Ciiq' - q^ + lC") + C^{q^) (4.16) 
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with /C3 = ci + C2 + C3 = Qi(0) -02(0) + Q3(0). 

If Ri satisfies the condition (|3.35p . the Birkhoffian (|4.16[) is dissipative. 
The function Eq^ : TM^ R and the dissipative 1-form satisfying (|2.15p are 
given by 

Xq,q) = J L^{e)q'dq^+ j L2{e)q'de + J C^iq' - q' + IC'){dq' - dq') + 
J C2{q')dq' + J C:i{q^)dq^ - j J C[{q' - q^ + lC')dq^dq^ (4.17) 

D = Ri{q^)dq^ (4.18) 

The equihbriuni points of the considered nonhnear network are the solutions of 
the system 

i?i(0) + Ci(gi -g2+/C=^) + C2(g^) = 

-Ci(gi -g2 + /c-^) + C3(g2) = (4.19) 

1) If i?i(0) = 0, and ii, L2, Ci, C2, C3 satisfies p.42p . then the equilibrium 
points are locally asymptotically stable. We define a Liapunov function V 

by 

V{qA)^E^^{q,q)~E^^{q,,0) (4.20) 

with Eq^ given by (|4.17p and q^ satisfying (|4.19p . The Hessian matrix of V at 
a equihbriuni point [q^. , 0) has the expression 



Hy(ge,0) = 



/ Li(0) \ 

£2(0) 

C[{ql,ql) + C!,{q',) -C[{qlql) 

\ -C[{qlql) C[{qlql) + C',{ql) I 

(4.21) 

Under the assumptions we made, this matrix is positive definite. The centers 
of the level curves of the function (|4.20p have the coordinates (ge , 0) , where q^ 
satisfies the system (|4.19p with i?i(0) = 0. Because the Birkhoffian (|4.16p is 
dissipative, it follows that dE^^ < 0, and therefore dV < 0. 
2) If i?i(0) ^ 0, but for aU x ^0 

X {Ri{x) - Ri{0)) > (4.22) 

and Li, L2, Ci, C2, C3 satisfies p.42p . then the equilibrium points are locally 
asymptotically stable. We define now a Liapunov function V by 

V(g,g) =EoJg,g)-EoJge,0) (4.23) 

with Eq^^ : TAIc R given by 

Eo^{q,q) = EQjq,q)+R,iO)q' (4.24) 
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Eo^iljQ) having the expression (|4.17p . The centers of the level curves of the 
function (|4.23[) have the coordinates (ge, 0), where qe satisfies the system (I4.19p . 
The Hessian matrix of the function V in (|4.23p has at the equilibrium point the 
same expression (|4.21l) . It remains to prove that dV < 0. This yields from the 
dissipativeness of the Birkhofhan (|4.16p . We consider the following dissipative 
1-form 

D ^[Ri{q') -R,{0)]dq' (4.25) 

In view of assumption (j4.22p . this vertical 1-form is indeed dissipative. One can 
easily check that for the function Eo^ (g, q) in (|4.24p and the dissipative 1-form 
in (|4.25p . the following identity is fulfilled 

2 2 

The Birkhoffian (|4.16p being dissipative we have dEiQ^ < 0, therefore dV < 0. 

Acknowledgement. This work was supported by IMAR through the con- 
tract of excellency CEx 06-11-12/ 25.07.06. 



dEo^ dEo^ 
—q H - 

dq^ dq^ 



g-' +Dj (g, q)q^ 



(4.26) 
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